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Abstract. We have observed an extended region sur- 
rounding the first discovered galactic superluminal source 
GRS 1915 + 105, seeking evidence of interaction between 
the relativistic ejecta of that object and the interstel- 
lar medium. We find two radio sources axisymmetri- 
cally aligned along the sub-arcsecond relativistic ejecta 
of GRS 1915 + 105 and roughly 17' distant from it, which 
coincide with the luminous IRAS sources 19124+1106 and 
19132+1035. We have observed these sources at centime- 
ter (VLA), millimeter (IRAM 30m), and infrared (ISO, 
UKIRT, ESO/MPI 2.2m) wavelengths in both line and 
continuum emission. At centimeter wavelengths a non- 
thermal jet-like feature aligned along the outflow axis is 
located adjacent to the inner edge of the southern source. 
Strong density enhancements are found in the millime- 
ter tracers CO and H 13 CO + at the positions of both 
sources and some of the morphology is reminiscent of 
shock-like interactions; however, linewidths are narrow. 
At infrared wavelengths strong hydrogen recombination 
lines and weak lines of molecular hydrogen are observed 
at the southern source. We discuss these results as possi- 
ble evidence of the sought-after interaction, both in terms 
of the regions undergoing ongoing shock-heating and in 
terms of them being locations of shock-induced star for- 
mation. The evidence for each of these is inconclusive. 
Millimeter line mapping of a portion of W 50 where the 
relativistic jets of the X-ray binary SS 433 interact with 
the interstellar medium shows roughly similar morphol- 
ogy as GRS 1915 + 105, suggesting that the phenomena 
observed at the IRAS sources may not be unusual for such 
a long distance interaction. 



Key words: Stars: individual: GRS 1915+105, SS 
433 - ISM: individual objects: IRAS 19124+1106, IRAS 
19132+1035 - ISM: jets and outflows - X-rays: stars 



1. Introduction 

The hard X-ray transient GRS 1915+105 was discovered 
in 1992 by the all-sky monitor WATCH on GRANAT 



(Castro-Tirado et ah, 1994). It has been intensively stud- 



ied since then in many wavebands, including the radio, 
where ejected plasma clouds with apparently superlu- 
minal velocities, a common extragalactic phenomenon, 
were found for t he first time in our Galaxy ( Mirabel and 
Rodriguez, 1991 ) . Material ejected from GRS 1915 + 105 
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would be expected to interact with the surrounding in- 
terstellar medium, providing an opportunity to study in 
detail for the first time such a relativistic interaction (for 
a review see Mirabel & Rodriguez (|1999|) ). 

In a program to look for these interactions we have de- 
tected two compact sources of bright radio emission, each 
coincident with a bright IRAS source near GRS 1915 + 
105. These sources are located axisymetrically with re- 
spect to GRS 1915 + 105 and at same position angle as its 
normal sub-arsec ejections, suggesting that they could be 
the zones of interaction between the ejecta and the inter- 
stellar medium as described in section ^. We have observed 
these two IRAS sources at near-infrared (ESO/MPI 2.2m, 
UKIRT), mid-infrared (ISO), millimeter (IRAM 30m) and 
centimeter (VLA) wavelengths. The observations are de- 
tailed in section |3| We discuss in section ^ the possi- 
bility of a physical association between the ejections of 
GRS 1915 + 105 and the two IRAS sources, comparing 
some of the results obtained in the putative interaction 
zones with new measurements of the interaction between 
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one of the jets of SS 443 and the interstellar medium. Some 
of the observations and results described here have been 
briefly described in Chat y (|1998| ), Chaty et al. ( ^000[ ) and 
Rodriguez and Mirabel (|1998|) (hereafter RM98). 



2. The context: two axisymmetric sources around 
GRS 1915+105 

In order to search for interactions involving the ener- 
getic and relativistic ejections of GRS 1915 + 105, the 
region surrounding GRS 1915 + 105 was observed at ra- 
dio wavelengths as described by RM98. This search was 
performed at A = 20 cm, using the Very Large Array 
(VLA) of NRAOQ, in its C-configuration, giving a resolu- 
tion of 15". The resulting map is shown in Figure 0. The 
region in the map referred to as G45.45+0.06 was misla- 
belled G45.46+0.06 in the original papers (e.g., RM98). 
Charac terist ics of this radio sourc e are given in Downes 
et al. (|1980D and in Feldt et al. (|l998| ). No evidence of 
jets or elongated clouds appears in the figure. However, 
two small radio continuum sources positioned nearly ax- 
isymmetrically with respect to GRS 1915+105 were found 
at angular separations of 17' each from GRS 1915 + 105 
(RM98). These small sources are coincident with the 
bright IRAS sources, 19124+1106 and 19132+1035. Their 
coordinates are given in Table [l]. The position angle of the 
line connecting the northwest source and GRS 1915 + 105 
is 157?9, and their separation is 16'6, the equivalent val- 
ues for the southeastern source and GRS 1915 + 105 are 
156P6 and 16'9. The position angles are very similar to 
the position angle of the sub-arcsec radio-ejections from 
GRS 1915+1 05 (~ 150°) (Mirabel & Rodriguez ( gggg ); 
Fender et al., (199?)). The angle between these ejections 
and the line of sight towards GRS 1915 + 105 is 70°, with 
the southeastern ejection approaching and the northwest- 



ern one receding (Mirabel and Rodriguez, 1994). 

The distance to GRS 1915+105 is crucial for determin- 
ing many of the physical parameters of the object (e.g., 
mass loss rate, velocity of ejecta, energetics). It wa s esti- 
mated to be 12.5 + 1 .5 kp c by Mirabel & Rodriguez ( 1994 ) 
and by Chaty et al. ( 1996 ). Recent VLB A obser vations are 
consistent with this value ( Dhawan et al., 2000|) . However, 
the uncertainty may be larger than 1.5 kpc, as the above 
lower limit was based on a distance to G45.45+0.06 of 



1 The National Radio Astronomy Observatory is operated 
by Associated Universities, Inc., under cooperative agreement 
with the USA National Science Foundation 



Source 


J2000.0 coord. 


gal. coord. 


GRS 1915+105 


a = 19 /l 15 m llf545 
8 = 10°56'44;'80 


l n =45? 40 
b n = -0?29 


IRAS 19124+1106 


a = 19 h 14 m 45f77 
8 = 11°12'06"4 


l n = 45? 54 
b u = -0?007 


IRAS 19132+1035 


a = 19 h 15 m 39fl3 
5 = 10°41'17"1 


I 11 = 45? 19 
b n = -0?44 



Table 1. Positions of GRS 1915+105, IRAS 
19124+1106 and IRAS 19132+1035. These coordi- 
nates are the positions of peak signal at 20 cm as observed 
by the VLA. 



GRS 1915+105 20 cm 1412.5 MHz 



11 30- 
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IRAS 19124+1 106 G 45. 8+0.13 



10 50 




8.8 kpc, whereas Feldt et al. ( 1998 ) have plac ed the H II 
region at only 6.6 kpc. Moreover, Fender et al. ( 1999 ), who 
derived an upper limit of 11.2 + 0.8 kpc, consider that the 
source is constrained to lie between 7 and 12 kpc. 

If the radio sources are at the same distances from the 
Sun as GRS 1915+105, and if we assume that the distance 
of this source is 12.5+1.5 kpc, their angular separations 
correspond to the distance of 60 pc from GRS 1915 + 105. 



1915 30 00 14 30 00 13 30 00 12 30 00 

RIGHT ASCENSION (B1950) 



Fig. 1. Map of the surroundings of GRS 1915 + 105, taken 
with the VLA-C at A = 20 cm. Contour levels are -3, 
3, 4, 6, 10, 15, 20, 30, 40, 60, 100, 200, 400 and 800 
x3 mjybeam -1 . The half power contour of the beam is 
shown in the bottom right corner. The arrows around 
GRS 1915 + 105 indicate the position angle of the sub- 
arc relativistic ejecta. 

3. Observations and results 

3.1. Centimeter 

High-resolution maps of the two radio sources have been 
obtained with the VLA. Details of these observations 
can be found in RM98. These maps, at three differ- 
ent wavelengths, 2, 6 and 20 cm, are shown in Fig- 
ure |^. The morphology of the northern source resem- 
bles that of a cometary HII region, but it also shows a 
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Source 


1.25 /im 


2.2 fim 


12 (im 


25 /im 


60 /im 


100 fim 


2 cm 


6 cm 


20 cm 




mag 


mag 


Jy 


Jy 


Jy 


Jy 


mjy 


mjy 


mjy 


IRAS 19124 + 1106 


17.9 ±0.1 


12.1 ± 0.1 


3.9 


19.6 


260.6 


581.5 


114 ±6 


130 ±4 


114 ±4 


IRAS 19132 + 1035 


17.45 ± 0.1 


10.7 ±0.1 


6.9 


34.0 


277.4 


488.8 


52 ±6 


63 ±4 


60 ±4 


IRAS 19132 + 1035 "jet" 














< 1 


2 


5 



Table 2. Flux densities of IRAS 19124+1106 and IRAS 19132+1035. 

All the flux densities with A from 12 /im to 20 cm come from RM98. We refer to Table || for the fluxes at mid-infrared 
wavelengths. 



bow shock-like structure to the South-East, e.g. towards 
GRS 1915 + 105. At the southern lobe a non-thermal jet 
extends to the northwest along the line between the source 
and GRS 1915 + 105. The flux densities of this jet are < 1, 
2 and 5 mjy at 2, 6 and 20 cm, respectively, indicating a 
spectral index of a = —0.8 (f v cx v a \ therefore noticeably 
different than that of the rest of the southern radio lobe, 
which exhibits thermal emission, as it can be seen in the 
Table ||. We also reported in Figure || a spectral index map 
of IRAS 19132+1035 made from the 20 and 6-cm maps. 
The southern lobe also shows a sharp edge to the south, 
which could be either a bow shock, or the ionization front 
of an H II region. 

Additional radio observations of these sources in the 
H92a recombination line and adjacent continuum at 
3.6 cm are detailed in RM98. The line strengths, profiles, 
and radial velocities are typical of H n regions at kine- 
matic distances of 7.4 ±1.4 kpc for the north lobe and 
6.0 ±1.4 kpc for the south lobe (distance errors are es- 
timated from the mean value for deviations of circular 



rotation of 12 kms~ given by Brand & Blitz ( 1993 ). We 
note that the distance uncertainties are large because of 
the low velocity resolution (±12 kms -1 ) of the H92 obser- 
vations (RM98). These estimates are most consiste nt wit h 
the distance to G45.45+0.06 given by Feldt et al. ([19980 . 
The radio luminosities of the sources derived using these 
distances are similar to those of H n regions powered by 
09.5 ZAMS and B0 ZAMS stars (RM98). If instead they 
are at the nominal distance of GRS 1915 + 105 hotter indi- 
vidual stars or two stars of the above spectral types would 
be required to power each source. However, the interpre- 
tation of these objects as normal Hn regions does not 
explain the presence of the non-thermal jet-like structure 
in IRAS 19132 + 1035. Furthermore, given the large un- 
certainty in the distance of GRS 1915 + 105, we can not 
a priori rule out the possibility that GRS 1915 + 105 is 
roughly at the kinematic distance of the H II regions. 

A search for OH maser emission was performed with 
the VLA in the vicinity of GRS 1915+105, at 1720 MHz. 
Such emission can be a signpost for interaction between a 
SNR and molecular gas (as in the case of IC443, e.g. De- 



noyer (1979)). No detection was obtained. However, this 



3.2. Infrared 
3.2.1. Near-infrared 

Broad-band imaging 

We imaged the IRAS sources in the J (1.25 fim) and 
K (2.2 fim) bands. The observations were made on 1997, 
April 5, with IRAC2b installed on the Max Planck In- 
stitute's 2.2 m telescope at the European Southern Ob- 
servatory (ESO). The IRAC2b camera, which contains a 
Rockwell 256x256 pixel Hg:Cd:Te NICMOS 3 array de- 
tector was mounted at the f/35 infrared adapter of the 
telescope. It was used with lens C, providing an image 
scale of 0.49 arcsec/pixel and a field of 136 x 136 arcsec 2 . 
The typical seeing for these observations was 1.2 arcsec. 
Each final image is the median of 5 frames, each exposed 
for 2 minutes. After taking each image of the object, an im- 
age of adjacent sky was taken, to allow subtraction of sky 
emission. The images were processed by removing bias and 
dark current, and applying a flat field correction. These 
steps were performed using IRAF procedures, in particu- 
lar the DAOPHOT package for the photometry in crowded 
fields. 

The images reveal a compact near-infrared counterpart 
at the position of IRAS 19124 + 1106, with J = 17.9 ±0.1 
mag and K = 12.1 ±0.1 mag, and a brighter and 
more extended counterpart of IRAS 19132 + 1035, with 
J = 17.45 ± 0.1 mag and K = 10.7 ± 0.1 mag. Figure 
^ shows the central portions of the images containing 
IRAS 19132+ 1035. 

Spectroscopy 

A low resolution K-band spectrum of IRAS 19132 + 
1035 was obtained with the facility instrument CGS4 on 
the United Kingdom Infrared 3.8 m Telescope (UKIRT) on 
UT 1997 July 13, as part of the UKIRT Service Program. 
The 75 1/mm grating was used in CGS4 together with a 
slit of width 1.2" to provide a spectral coverage of 0.67 pan. 
at a resolution of 0.0026 /im (370 kms -1 at 2.15 fim) on 
the 256 x 256 array of InSb detectors. The slit was oriented 



does not imply necessarily that there is no interaction (Ro- 
driguez, Goss, Mirabel, private communication). 



2 Based on observations collected at the European Southern 
Observatory, Chile. 
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Fig. 2. Maps of the two continuum radio sources 
IRAS 19132+ 1035 (left) and IRAS 19124+ 1106 (right), 
acquired with the VLA respectively in the D-configuration 
for the A = 2 cm map (top), in the C-configuration for 
the A = 6 cm map (middle) and in the B-configuration 
for the A = 20 cm map (bottom). Contour levels are 
-4, 4, 6, 8, 10, 15, 20, 40, 60, 100, 200, 300 and 400x the rms 
noise of 0.05 mJybeam -1 . The half power contour of the 
beam, with diameter of 4", is shown in each bottom left 
corner. 




19 h i5"39?a0 3^30 3^20 3<)=10 3^00 39!<J0 38?8Q 38*70 



10"41 '10' 
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Fig. 4. J (top) and K (bottom) band images of the source 
IRAS 19132 + 1035, taken with IRAC2b on the 2.2 m of 
the ESO. The coordinates are for the J2000 equinox. The 
ADU flux contour levels are in a logarithmic scale, each 
separated by a factor \/2, the first one is equal respectively 
to 0.864 for the J band and to 25.501 for the K band. 



at a position angle of 148° and the telescope was moved 
to center the peak signal from the infrared counterpart in 
one row of the array. Observations were obtained in the 
standard stare / nod-along-slit mode. The total integra- 
tion time was 8 minutes. A near simultaneous spectrum of 
the F3V star HR 6987 (T=6700 K, K=4.50 assumed) was 
obtained at the same airmass; in order to remove telluric 
absorption lines (the Brackett 7 (7-4) absorption line, with 
a central depth of 0.88 of the continuum at this spectral 
resolution, was artificially removed from the spectrum of 
the comparison star prior to ratioing). Wavelength cali- 



bration was obtained from observations of an argon arc 
lamp and is accurate to 0.0005 ^m. 

The spectrum of IRAS 19132 + 1035 is shown in Fig- 
ure ^. It is the sum of the spectra of three adjacent rows, 
covering an area of 1.2 x 3.6 arcsec. Most of the flux is con- 
tained in the central row. The spectrum exhibits a very red 
continuum and a number of emission lines. Most promi- 
nent are three recombination lines of atomic hydrogen: a 
very strong Paschen a (4-3) line at 1.876 /mi, and weaker 
Brackett 7 (7-4) and Brackett S (8-4) lines, at 2.166 fim 
and 1.945 /im, respectively. The Pa a and Br S lines are 
uncorrected for the absorption due to the same line in the 
comparison star and both it and the Br S line occur near 
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3.2.2. Mid-infrared 




Both infrared sources were imaged in six wavebands from 
7 to 15 fim with the infrared camera ISOCAM^ QCcsarsky| 



et al., 1996), on the Infrared Space Observatory (ISO) 



Wavelength (fJ.ni) 

Fig. 5. K-band spectrum of IRAS 19132+ 1035 obtained 
with CGS 4 on UKIRT. The upper trace is multiplied by 
a factor of 3.33 to allow an easier reading of the weak- 
est lines. Identifications and positions of detected and 
marginal lines are indicated. 



strong telluric absorption lines of water vapor; thus their 
strengths and profiles are subject to systematic errors. 

The peak of the Br 7 emission line occurs at 
+45 ± 70 kms -1 in the Local Standard of Rest (LSR). 
The Full Width at Half Maximum (FWHM) of the pro- 
file is only marginally broader than the instrumental res- 
olution, but the profile appears to exhibit a weak high- 
velocity wing, displaced towards the red. No such wing 
is present in the profile of either of the other two atomic 
hydrogen lines. We note that lines displaced towards the 
blue are much more common in the H II regions, because 
the front part of an expanding H II region is more easily 
detected. 

Two fainter emission lines also are visible: the sin- 
glet He I 2P — 2S transition at 2.059 /im and the 
H 2 1 - 5(1) line at 2.122 fjtm. The flux in the H 2 line 
is ~ 1 x 10~ 17 W m~ 2 . There is also evidence for the 
presence of the H2 2 — 1 S(l) line at 2.224 /im, at roughly 
one- fourth the strength of the 1-0 line, which suggests that 
most of the excitation of Hi is due to absorbed UV ra- 
diation rather than collisions. The 1-0 line appears partly 
resolved with a deconvolved FWHM of ~ 500 kms -1 , 
but new measurements are required to confirm this. The 
signal-to-noise ratio of the 2-1 line is too low to determine 
if its profile is resolved. 



satellite mission as part of the Guaranteed and Open Time 
Programme. The observations were performed on 1996 
April 28 and on 1997 October 20. The log of the observa- 
tions and the derived fluxes in all bands are found in Table 
||. The images were taken through the Large Width (LW) 
filters at angular resolutions given in this Table. Data re- 
duction used the CIA package and included subtraction 
of the dark current, suppression of the cosmic ray impacts 
by a multiresolution median method, correction of the de- 
tector transient behaviour, flat field correction and when 
necessary a distortion correction (see Starck et al. ( 1999 ) 
and references therein). 

The observation at 7 fim of the south lobe, shown at 
the bottom of Figure |^, is particularly striking, because 
it shows thermal emission, presumably from heated dust, 
exactly coincident with the radio counterpart. One compo- 
nent of this emission is localized, between the maximum 
and the jet feature seen at radio wavelengths, on a line 
connecting GRS 1915+105 and IRAS 19132+ 1035. The 
other component is more extended, and closely resembles 
the remainder of the radio image with its bow-shock mor- 
phology at the south-east edge. The 15 /im image of the 
north lobe is shown at the top of Figure ||. The morphol- 
ogy is similar to the radio images, resembling a common 
cometary H II region, but with two bow shock-like struc- 
tures facing each other. 

3.3. Millimeter 



We used the 30 m radio telescope of the Instituto de Ra- 
dio Astronomfa Milimetrica (IRAM0) to observe molecu- 
lar lines which are good density tracers and are diagnostics 
of shock excitation and chemistry. 

The observations, totalling 70 hours, took place during 
1997 April 17 - 24, November 29 - December 3 , and 
during 1997 December 31 - 1998 January 5. The observed 
molecules and transitions, together with their respective 
frequencies were: 12 CO (J = 2 - 1) at u = 230.5 GHz, 
13 CO (J = 2 - 1) at v = 220.4 GHz, H 13 CO+(J = 1-0) 
at v = 89.2 GHz, SiO (J = 2 - 1, v = 0) at v = 86.8 
GHz, SiO (J = 3 - 2,v = 0) at v = 130.3 GHz, SiO (J = 



3 ISOCAM was constructed under the scientific direction and 
the technical expertise of the Service d'Astrophysique of the 
CEA/Saclay. The ISOCAM Consortium is led by the PI C. 
Cesarsky. 

4 IRAM is an European institute for millimeter astronomy, 
founded by the French Centre National de la Recherche Sci- 
cntifique (CNRS) and the German Max-Planck-Gesellschaft 
(MPG), in collaboration with the Spanish Instituto Geografico 
Nacional (IGN). 
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Source 


Date 


Filter 


PFOV 


Flux 


IRAS 19132 + 1035 


28/04/96 


LW2 (5 - 


8.5 (im) 


3" 


13.2 ±2.6 Jy 


IRAS 19132 + 1035 


28/04/96 


LW7 (8.7 - 


10.7 /im) 


6" 


5.9 ± 1.2 Jy 


IRAS 19132 + 1035 


28/04/96 


LW3 (12 - 


■ 18 fim) 


6" 


10.7 ± 2.2 Jy 


IRAS 19132 + 1035 


20/10/97 


LW10 (8 - 


- 15 /im) 


1.5" 


6.3 ± 1.2 Jy 


IRAS 19124 + 1106 


20/10/97 


LW9 (14- 


- 16 /-im) 


3" 


3.3 ± 0.66 Jy 



Table 3. Log, characteristics and fluxes of the sources from ISO observations. 

PFOV = Pixel Field Of View. 



5 - 4,<u = 0) at v = 217.1 GHz and CS (J = 2 - 1) at 
v = 98.0 GHz. Details are given in Table f|. 

We used position switching during the observations, 
with the off-source position located (—500",— 1200") from 
the position of GRS 1915 + 105. Tests demonstrated that 
the off position contained no significant line emission. The 
treatment and analysis of the millimetric data were stan- 
dard, consisting for each spectrum of a subtraction of a 
polynomial baseline of first or second order, in order to 
remove instrumental fringes. 

3.3.1. IRAS 19132+ 1035 

For each of the 12 CO, 13 CO, H 13 CO + and SiO 2-1 
transitions, we obtained 23 spectra at positions separated 
by 6", along the line joining GRS 1915 + 105 to the 
South lobe and crossing the maximum of the centimetric 
continuum shown at the bottom left panel of Figure 0; 
this maximum was defined as the (0,0) position. Integra- 
tion times were 15 minutes per point, except at the (0,0) 
position where it was 45 minutes, (-13, 20) where it was 
135 minutes, and (-10, 15) where it was 90 minutes. The 
two latter locations correspond to the jet-like feature. In 
addition we obtained spectra at eight locations adjacent 
to the jet. The observations of CS, SiO 3-2 and SiO 
5-4 each consisted of 15 spectra at 6" separation along 
the same line, with integration times of 15 minutes per 
position with the exception of (0,0) and (10,-15) where 
the times were 30 and 7 minutes, respectively. 

Results for 12 CO, 13 CO, H 13 CO + and CS 

Figure [t] contains position-velocity plots with an 
expanded velocity scale for the observed lines of 12 CO, 
13 CO, H 13 CO + and CS. For all transitions the line peak 
occurs at an LSR velocity of 67 kms - , close to the 
velocity of 75.7 kms -1 of the H92a recombination line 
observations as seen by RM98 with much lower velocity 
resolution. The intensity distributions of the millimeter 
lines are not symmetric, each showing a sharper edge to 
the south (positive RA offset), as seen in the VLA images, 
and a slightly shifted peak velocity. The CS transition, 
which is a high density tracer, exhibits a peak shifted to 
the northwest (more negative RA offsets), towards the 
jet-like feature. 




Fig. 8. Spectrum of the SiO 2-1 line at IRAS 19132+1035, 
obtained at the offset (a, 5) = (-10", +14"). In this and the 
following spectra, the y-axis is the antenna temperature 
in K and the x-axis is the radial velocity with respect to 
the LSR. 

SiO 

The SiO 2-1 transition was detected only close to 
the position (-10, 14); the spectrum at that location is 
shown in Figure ||. The line is weak, but its velocity, 
67.6 ± 0.5 kms -1 and width, 4.4 ± 1.2 kms -1 are similar 
to those seen for much stronger lines and give us confi- 
dence that it is real. The 3-2 transition was not detected 
at any individual position, but the average of the spectra 
between (0,0) and (-17,25) yields a 3 a detection of a nar- 
row (1.0+0.3 kms -1 ) line at 67.3 kms -1 ±0.2, exhibiting 
an antenna temperature of T\ = 0.044 K and an intensity 
of 0.05 + 0.01 K km s -1 (Figure g). The SiO 5-4 transition 
was not detected at any individual location or when all of 
the spectra were averaged. 

Discussion 

Figure ^ demonstrates clearly that brightness maxima 
for transitions with higher critical densities are displaced 
increasingly towards the northwest where the non-thermal 
jet-like structure is located. What little SiO emission is de- 
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Abbr 


Transition 






Frequency 


Receiver 


Banks 


Beam 


12 CO 


12 CO (J = 2 


-1) 




v = 230.537 990 GHz 


230g2 


b44 b34 


10.9" 


13 co 


13 CO (J = 2 


-1) 




v = 220.398 686 GHz 


230gl 


b43 b33 b20 


11.4" 


H 13 CO + 


H 13 CO+(J = 


1 - 


0) 


v = 89.188 523 GHz 


3mml 


b42 b32 


28.2" 


CS 


CS (J = 2 - 


1) 




v = 97.980 968 GHz 


3mml 


b41 b31 


25.7" 


SiO 2-1 


SiO (J = 2 - 


l,v 


= 0) 


v = 86.846 891 GHz 


3mm2 


b41 b31 


29" 


SiO 3-2 


SiO (J = 3 - 


2,v 


= 0) 


v = 130.268 702 GHz 


2mm 


b42 b32 


19.4" 


SiO 5-4 


SiO (J = 5 - 


4,v 


= 0) 


v = 217.104 935 GHz 


230gl 


b43 b33 b20 


11.6" 



Table 4. Frequencies of the transitions and receivers used at IRAM. The filter bank b3*, divided into sections 
b31, b32, b33 and b34, provides a low resolution of 1 MHz, and contains 4 x 256 channels, for a total bandwidth of 1 
GHz; the fhterbank b20 provides a high resolution of 100 kHz, and contains 2 x 128 channels, for a total bandwidth of 
25.3 MHz; the autocorrelator b4*, divided into sections b41, b42, b43 and b44 offers resolutions ranging from 10 kHz 
to 1.25 MHz, a variable number of channels, and a variable total bandwidth from 20 MHz to 1 GHz. 




50 100 

Velocity (km/s) 



Fig. 9. Spectrum of the SiO 3-2 line at IRAS 19132+1035, 
averaged over all positions from (0,0) to (-17;25). 



tectcd also originates largely in this region. The very dens- 
est region could be the main interaction zone between the 
putative jet and the ambient medium. However, the nar- 
row lines that are observed are difficult to explain in the 
above scenario, as an energetic shock might be expected 
to lead to line emission over a much wider range of veloc- 
ities than is observed. If a shock is present, it apparently 
must be of very low velocity. 

3.3.2. IRAS 19124+1106 

For the 12 CO, 13 CO, H 13 CO + and SiO 2-1 transitions, we 
acquired 24 spectra at positions separated by 6", along 
the line joining GRS 1915 + 105 to the northern lobe. As 
in the case of the southern lobe, the reference position 
(0,0) is defined as the position of maximum radio contin- 
uum signal in the northern lobe, as observed by the VLA. 
The integration time at each position was 5 minutes, 



except at the reference position, where it was 10 minutes. 
Each of the above transitions was detected over a wide 
range of positions. Spectra of the CS 2-1, SiO 3-2 and SiO 
5-4 transitions were obtained only at the (0,0) position, 
each with an integration time of 15 minutes. None of the 
three SiO transitions were detected at any position at an 
upper limit of 0.05 K, or when all positions were averaged. 

Results for 12 CO, 13 CO and H 13 CO + 

Position velocity diagrams for 12 CO, 13 CO and 
H 13 CO + are presented using an expanded velocity scale 
in Figure Two strong velocity components are 

readily apparent in all three transitions. For 12 CO at the 
(0,0) position, the strongest component is centered at 
53.95 ± 0.03 kms" 1 , with a width of 4.05 ± 0.09 kms -1 , 
an intensity of 12.9 K, and an integrated brightness 
of 55.5 ± 0.09 Kkms -1 . For 13 CO, this component 
is peaked at 55.10 ± 0.04 kms - , and has a width 
of 3.88 ± 0.10 kms -1 , an intensity of 7.32 K, and a 
integrated brightness of 30.3 ± 0.60 Kkms -1 . A second 
velocity component occurs at 59.5 ± 0.1 kms -1 with a 
width of 3.9 ± 0.2 kms -1 , an intensity of 3.11 K and an 
integrated brightness of 12.8 ± 0.6 Kkms -1 at the 12 CO 
(0,0) position. The velocities of these two components 
are close to that of the H92a recombination line peak 
at 57.3 kms -1 (RM98). A third velocity component, at 
~ 6 kms -1 is visible in the northwestern part of this 
cloud for CO and in the southeastern part for H 13 CO+. 
Given the very different velocity the line emission proba- 
bly originates in a foreground or background cloud. 

Results for CS 

The single spectrum obtained at (0,0), shown in 
Figure O, contains a strong component centered at 
56.17 ± 0.03 kms -1 , of width 3.039 ± 0.07 kms -1 , 
intensity 0.937 K and integrated brightness 
3.03 ± 0.06 Kkms -1 . Adjacent to it is a weaker 
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ered by one or more hot stars, and (3) the cometary mor- 
phology of the northern source, IRAS 19124+1106, which 
is common in star forming regions. 

Although the physical alignment of the IRAS sources 
could be a background coincidence, the probability of this 
is low in view of the small number of IRAS sources, the 
existence of a jet at IRAS 19132 + 1035, and the loca- 
tions of these sources within the Galaxy. It remains possi- 
ble that the jet at IRAS 19132 + 1035 is an extragalactic 
background source. This should be tested by radio obser- 
vations at higher resolution. Because an association be- 
tween GRS 1915 + 105 and these objects may be real we 
explore in the following the potential characteristics of 
interactions of ejecta from GRS 1915 + 105 with the inter- 
stellar gas and compare these examples with observations 
of another somewhat similar source. 



Fig. 11. IRAS 19124+1106, CS, b41, spectrum taken at 
the position (0,0). 

component at 59.91 ± 0.11 kms -1 , with width 
1.589 ± 0.22 kms -1 , intensity 0.298 K and integrated 
brightness 0.504 ± 0.07 Kkms -1 . These velocities are 
similar to those detected in other molecules and in H92o; 
(RM98). 

Discussion 

As already noted by RM98 from the centimeter obser- 
vations, the morphology of the northern source is that of a 
cometary H II region. The millimeter wave spectroscopy, 
which reveals two velocity components, is consistent with 
this morphology. 

4. Discussion 

The new data do not definitively prove or disprove an as- 
sociation between GRS 1915+105 and either of the IRAS 
sources. The strongest lines of evidence supporting an as- 
sociation are (1) the axisymmetric locations of the two 
sources at nearly the same position angle as the recent 
sub-arcsecond ejections observed at radio wavelengths, (2) 
the spatial coincidence of the non-thermal radio jet with 
the inner edge of the southern source, IRAS 19132 + 1035, 
as well as the orientation of this jet along the axis, (3) the 
location of the highest densities in IRAS 19124+1106 and 
IRAS 19132 + 1035 on the sides closest to GRS 1915 + 105 
and, in the case of IRAS 19132 + 1035, close to the non- 
thermal jet, and (4) the bow shock - like structure in 
the portion of IRAS 19132 + 1035 most distant from 
GRS 1915 + 105. Some of these lines of evidence are asso- 
ciated only with the southern source. Evidence against an 
association includes (1) the lack of detected high velocity 
gas at the IRAS sources, (2) the luminosities of the IRAS 
sources, which are consistent with each of them being pow- 



4-1- Production of a large scale cavity 

The angular separation between each IRAS sources and 
GRS 1915 + 105 is 17'. The measured proper motions of 
the most recent ejections from GRS 1915 + 105 are respec- 
tively for the approaching and receding ejecta in the inter- 
val Happ — [17.6, 23.6] ±0.5 mas/d a nd \i ve c = [9.0, 10.0] ± 
0.5 m as/d (Mirabel & Rodriguez ( |l994| ), Fender et al. 
( 1999| ) ) . This implies typical travel times to reach the lobes 
of 140 years for the approaching ejecta and of 300 years for 
the receding ejecta. It is therefore interesting to note that 
the impact due to the South lobe with the southern IRAS 
source would appear to occur much earlier than that of the 
North lobe with the northern IRAS source. Evidence for 
precession of the jets during the last five years is marginal 
(Rodriguez and Mirabel, 1999), with the maximum angle 
no more than 10 degrees. Thus a close alignment of cur- 
rently observed jets close to GRS 1915 + 105 with those 
from long past outbursts that may be responsible for the 
IRAS sources may not be a coincidence. 

It is perhaps surprising that the separations from 
GRS 1915 + 105 of the two possible interaction zones are 
identical to better than two percent. These are determined 
not only by properties of the outflow from GRS 1915+105, 
but also by the density and uniformity of the interstel- 
lar medium. In order for ejected matter to travel up to 
60 pc, a large scale cavity must have been created around 
GRS 1915 + 105. Over time this cavity could tend to ex- 
tend an equal distance in each direction. The probable 
existence of a cavity around GRS 1915 + 105 already has 
been pointed out by Mirabel and Rodriguez ( |i996l ). The 
plasma clouds ejected during the last few years exhibit bal- 
listic motions (i.e., with no evidence of deceleration), indi- 
cating that the density in the interstellar medium through 
which they have passed is less than 0.05 protons cm~ 3 
(Mirabel & Rodriguez, 1996). These clouds have been ob- 
served until their emission faded out or became unobserv- 
able with the available VLA configuration, typically at 
distances of one arcsecond from the central source. 
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Although direct evidence for a cavity is limited to a 
region very close to GRS 1915 + 105, it is likely that the 
cavity extends well beyond one arcsecond from it. Apart 
from the two IRAS sources, no phenomena that might be 
associated with deceleration of recent ejecta have been ob- 
served from 1" angular separation from GRS 1915 + 105 
to 17'. In order to produce a cavity with a linear dimen- 
sion of many tens of parsecs, even a relativistic source 
must be active for a much longer period of time than 
GRS 1915 + 105 has been observ ed. From EXOSAT mea- 
surements ( Reynolds et al., 1999 ) GRS 1915 + 105 was ac- 



tive in 1985 but much fainter than it has been since 1992. 
For the second discovered superluminal galactic source, 
GRO J1655 - 40, Kolb et al. ( |1997| ) suggested that the 
activity of the source is caused by the companion crossing 
the Hertzsprung gap, a transient phase lasting ~1 Myr. 
The companion of GRS 1915 + 105 ( [Marti et al., 2000D 
is believed to be more massive than the companion of 
GRO J1655 - 40 ( jBhahbaz et al., 1999[ ), and thus the du- 
ration of the active phase would be much shorter, although 
still sufficiently long, especially since GRS 1915 + 105 has 
been much more active than GRO J1655 — 40. 



4-2. Interaction between ejecta and ambient material 

Considerable heating of the gas must occur at the ends of 
the cavities, where the ejecta, whether or not partially de- 
celerated en route, impact ambient material. One would 
expect that much of heating would occur via shock ex- 
citation. For the case of GRS 1915 + 105 the speed and 
density of the incident material and the relative amounts 
of atomic and molecular gas (both in the incident and im- 
pacted material) are not known, and thus the spectrum 
of the shock-heated gas is uncertain. The southern IRAS 
source shows strong infrared lines of atomic hydrogen and 
weak infrared line emission from molecular hydrogen, very 
unlike the emission observed in molecular clouds impacted 
by protostellar winds having speeds of up to 100 kms -1 , 
in which lines of molecular hydrogen dominate and hy- 
drogen recombination lines are nearly non-existent (e.g. 
Geballe & Garden Q1987D ). 

We note that if the non-thermal feature near the 
southern lobe is an element of the interaction involving 
the ejecta of GRS 1915 + 105, the collimation of the jet 



is U 



length j e t 



w ^ 100, ten times larger than 



width j 

the most highly collimated Herbig-Haro objects (see e.g. 
Bachiller ( |1996| )). This implies a jet opening angle of 
< 0.5° which is consistent with the limit of < 8° derived 



by Fender et al. (1999) based on observations of the core, 
but is more highly constraining. This jet opening angle 
seems to be in agreement with the jet full opening angle 



as function of distance from the core for M87 (Junor 
et al., 1999[ ), if 

we take into account the scaling factor 



between galactic and extragalactic black holes. Therefore, 
as for extragalactic jets, it would require in the galactic 
case ongoing confinement of the jet at large distances 



from the source, ruling out free expansion of a relativistic 
gas. 

Induced star formation ? 

The similarity of the luminosities of the IRAS 
sources to those of compact H II regions containing 
one or more massive stars, as well as the spectrum of 
IRAS 19132 + 1035, which shows strong hydrogen recom- 
bination line emission, suggest the possibility that ejecta 
from GRS 1915 + 105 have induced massive star forma- 
tion at the locations of the IRAS sources, via compression 
of the interstellar gas. In this scenario the non-thermal 
radio jet observed at IRAS 19132 + 1035 might be in- 
terpreted as a Herbig-Haro-like feature, the result of a 
protostellar wind breaking out of the natal cloud. How- 
ever, given the length of time for star formation to pro- 
ceed to this phase, such an explanation would require ac- 
tivity from GRS 1915 + 105 at a much earlier time than 
mentioned above. Moreover, there appears to be no rea- 
son for star formation at IRAS 19132 + 1035 to produce a 
close coincidence between the position angle of the radio 
jet at IRAS 19132 + 1035 and that of the line connect- 
ing GRS 1915 + 105 and IRAS 19132 + 1035, as is ob- 
served. Finally, the radio emission from Herbig-Haro jets 



usually is thermal (free-free; Rodriguez (1999)). In view 
of these arguments, we do not regard induced star for- 
mation to be a likely explanation for the luminosities of 
IRAS 19124+ 1106 and IRAS 19132 + 1035. Dubner et al. 
(1998) note that in the case of the relativistic jets from 
SS 433 there is no evidence of induced star formation in 
the impacted gas. 

4-3. Comparison with SS 433 

In view of the inconclusive evidence linking the IRAS 
sources with activity originating at GRS 1915 + 105, it is 
perhaps useful to examine other possible examples of such 
an interaction. Like GRS 1915 + 105, the famous X-ray bi- 
nary SS 433 ejects beams of material at relativistic speeds. 
This object is located inside the radio shell/supernova 
remnant W50, which exhibits two lateral extensions with 
dimensions of tens of parsecs (Dubner et al., 1998). The 



morphology can be attributed to continuous ejection of 
magnetic field and high-energy particles from the central 
source. Sub-arcsecond jets are present at SS 433, with 5 
orders of magnitude difference between their extents and 
those of the far lobes. The power injected int o the se jets 
is ~ 10 39 ergs -1 according to Dubner et al. (1998 ) , who 
also estimate that the kinetic energy transferred to the 
surroundings of SS 433 during the last ~ 2 x 10 4 years 
amounts to 2 x 10 51 erg. 

Major ejection events in GRS 1915+105 are more 
sporadic than in SS 433, and therefore the average ki- 
netic energy injected into the surroundings over long pe- 
riods of time may be smaller (Mirabel and Rodriguez 
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1999). The kinetic energy of the 1994 March 19 event in 
GRS 1915 + 105 was - 10 43 erg ( [Rodriguez and Mirabel" 
1999()7]GRS 1915 + 105 would need to exhibit one ejection 
event per hour similar to this one in order to match SS 443. 



A more likely rate is one event per month (Rodriguez and 
Mirabel, 1999), in which case ~ 2 x 10 7 years would be 
required to equal the energy output of SS 443 in 2 x 10 4 
years. However, we note also that the jets emanating from 
GRS 1915 + 105 alwa ys are compact close to the source 
which indicates that a continuous 
105 is occuring. In 



(phawan et al., 20001) , 
injection of energy from GRS 1915 
addition to the large events described above, there are 
also smaller ones, observed at X-ray, radio and infrared 
wavelengths, where the mechanical luminosities are in the 
range 10 37 — 10 39 erg s -1 , and where the synchrotron emis- 
sion isseen up to the infrared wave lengths (e.g. Mirabel et 
al. (tL998| ), Fender & Pooley fl2000| )). A comparison of the 
energetics of GRS 1915 + 105 with SS 433 and a Herbig- 
Haro object is reported in Table ||. 

Despite the differences between GRS 1915 + 105 and 
SS 433, each may have created a cavity of similar linear 
extent. In view of the uncertainty as to what effects might 
be observable at the putative ends of the cavity produced 
by GRS 1915 + 105, we have observed the western tip of 
W50 where the ejecta of SS 443 and the shell remnant in- 
teract, in some of the same millimeter lines of 12 CO, 13 CO, 
H 13 CO+ and SiO as were observed for the IRAS sources. 
The observed region in W50 is at the constant declina- 
tion 05°00'00" and right ascensions between 19 h 05 m 27 s 
and 19' l 06 m 24 s (B1950) (see Dubner et al. ( ]199S| ), fig- 
ures la and lb). The off position was 15 arcminutes east 
and north of the reference point located at 19 /l 06 m 00 ;! ; 
05°00'00". Measurements at W50 were separated by 24", 
with integration times of 5 minutes each, except at (0,0) 
where the integration time was 10 minutes. 

The results are shown in Figure \L2 over the relevant 
velocity range. Emission from 12 CO, ^CO and H 13 CO + 
are found roughly 200 - 300" west of the reference po- 
sition, at the westernmost edge of the radio shell (Dub- 
ner et al. ( |1998| ), Figures la and lb). The emission from 
the high density tracer H 13 CO + is the most compact and 
probably indicates the location of the strongest interac- 
tion. The line profiles of all three species are asymmetric 
but narrow. SiO 2-1 emission was not detected. These 
results are very similar to those for IRAS 19132 + 1035. 
These similarities between SS 433/W50 and GRS 1915 + 
105/IRAS 19132 + 1035 and IRAS 19124 + 1106 suggest 
that high density, low velocity molecular clumps may not 
be unusual products of the interaction of the ejecta of a 
distant energetic source and the surrounding medium. 



5. Conclusions 

We have performed extensive multi-wavelength observa- 
tions from IR to radio of two radio/IRAS sources axisy- 
metrically located with respect to GRS 1915 + 105 and 





HH 


SS 433 


GRS 1915 


Velocity 


1 - 500 km/s 


0.26 c 


0.92 c 


Separation 


lau — lOpc 


50 pc 


60 pc 


Range 


10 6 


10 5 


10 4 - 10 5 


Collimation 


2-10 


> 12 


~ 100? 


Kinetic Pow. 


10 32 ergs^ 1 


10 39 ergs" 1 


10 36 - 39 ergs" 1 



Table 5. Characteristics of jets from different 
sources. HH = Herbig-Haro. The range is the ratio of 
distance between the far lobes and the subarcsec jets, the 
collimation factor is the ratio between the length of the 
jet and its width, and the kinetic power of the ejections is 
a time average. 



aligned with the position angle of the subarcsec jets. 
The northern source has the morphology of a common 
cometary H II region. The observations of the southern 
cloud reveal a collimated non-thermal structure which 
may be the far end of a jet of material from GRS 1915 + 
105. Both sources contain dense clumps of molecular ma- 
terial. Overall the evidence for these regions being inter- 
action zones seems inconclusive. An abbreviated study of 
what may be a very similar interaction between the ejecta 
of SS 433 and the surrounding interstellar medium reveals 
strikingly similar (and also apparently inconclusive) phe- 
nomena, perhaps indicating that such phenomena might 
be produced by sources which via their violent behavior 
have created large cavities over long periods of time. 
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Fig. 3. Spectral index map of IRAS 19132+1035 made 
from the 20 and 6-cm maps. The color coding for the spec- 
tral index is given at the bottom of the figure. Note the 
negative, non-thermal spectral index of the jet feature to 
the northwest. 
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Fig. 6. Top: ISO map of the source IRAS 19124 + 1106, 
at A = 15 /im, taken with the LW9 filter. Superimposed 
are the 20cm radio contours at 0.4, 0.6, 1, 2, 5, 10 and 15 
mJy. Bottom: ISO map of the source IRAS 19132 + 1035, 
at A = 7 /Ltm, taken with the LW2 filter. Superimposed are 
the 20cm radio contours at the levels 0.2, 0.4, 0.7, 1.2, 1.8, 
2.5, 3, 4 and 5 mJy. 
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Fig. 7. Observations of IRAS 19132+1035. Offsets are 
relative to the position of maximum radio emission ob- 
served at the VLA. The black contours are antennae iso- 
temperature: 12 CO: = — 1, and from 1 to 2QK sepa- 
rated by an interval of IK; 13 CO: = —1, and from 1 to 
11 K separated by an interval of 1 K; H 13 CO+: = -1, 
and from 0.2 to 2.1K separated by an interval of 0.1K- 
; CS: = — 1, and from 0.2 to 2.2K separated by an 
interval of 0.2i^. 



Fig. 10. Observations of IRAS 19124+1106. Offsets are 
relative to the position of maximum radio emission ob- 
served at the VLA. The black contours are antennae iso- 
temperatures: 12 CO: T^ = -1, and from 1 to 14 K sepa- 
rated by an interval of 1 K; 13 CO: -1 and from 0.5 to 8 K 
separated by an interval of 0.5 K; H 13 CO + : -1 0.3 0.4 and 
from 0.5 to 0.8 K separated by an interval of 0.05 K. 
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Fig. 12. Observations of W50, the supernova remnant 
shell surrounding SS 433. Offsets are relative to the posi- 
tion of the central source. The black contours are antennae 
iso-temperatures: 12 CO: = -1 and from 1 to 12 K sep- 
arated by 1 K; 13 CO: -1 and from 0.5 to 5 K separated by 
1 K; H 13 CO+: -1 and from 0.2 to 1 K separated by 0.1 K. 



